Dynamical analysis of a damped harmonic forced duffing oscillator with time delay

This paper is concerned with a time-delayed controller of a damped nonlinear excited Duffing oscillator (DO). Since time-delayed techniques have recently been the focus of numerous studies, the topic of this investigation is quite contemporary. Therefore, time delays of position and velocity are utilized to reduce the nonlinear oscillation of the model under consideration. A much supplementary precise approximate solution is achieved using an advanced Homotopy perturbation method (HPM). The temporal variation of this solution is graphed for different amounts of the employed factors. The organization of the model is verified through a comparison between the plots of the estimated solution and the numerical one which is obtained utilizing the fourth order Runge–Kutta technique (RK4). The outcomes show that the improved HPM is appropriate for a variety of damped nonlinear oscillators since it minimizes the error of the solution while increasing the validation variety. Furthermore, it presents a potential model that deals with a diversity of nonlinear problems. The multiple scales homotopy technique is used to achieve an estimated formula for the suggested time-delayed structure. The controlling nonlinear algebraic equation for the amplitude oscillation at the steady state is gained. The effectiveness of the proposed controller, the time delays impact, controller gains, and feedback gains have been investigated. The realized outcomes show that the controller performance is influenced by the total of the product of the control and feedback gains, in addition to the time delays in the control loop. The analytical and numerical calculations reveal that for certain amounts of the control and feedback signal improvement, the suggested controller could completely reduce the system vibrations. The obtained outcomes are considered novel, in which the used methods are applied on the DO with time-delay. The increase of the time delay parameter leads to a stable case for the DO, which is in harmony with the influence of this parameter. This drawn curves show that the system reaches a stable fixed point which assert the presented discussion.

. Dot signifies the derivative of y with regard to t. ω Root of the natural frequency of the scheme. µ Pure linear damping coefficient. α Cubic nonlinear Duffing coefficient, α > 0 , and α < 0 are hardening and softening spring, respectively. β Quintic nonlinear Duffing coefficient. F Force's amplitude.
Force's frequency. www.nature.com/scientificreports/ In reality, oscillation is a regular occurrence. It has the potential to both cause damage and yield benefits in life. Whenever the oscillation reaches a particular threshold, it can cause damage to human and animal health, keeping the human body discomfort and preventing technology from working regularly. Vibrations have advantages where they can be transformed into other forms such as types of energy, create lovely music, send out signals, and so on. This will have a significant impact on the advancement of human living standards. Consequently, vibration controllers and employment have been developed as fundamental technology and challenge in the mechanical development and production of current high-tech machinery. The relevance of time-delayed differential equations has led to the development of a substantial area of nonlinear analysis with extensive applications in biology, economics, physics, nuclear reactor, chemistry, and natural science fields. In powerful vibrating controllers, the time delay is necessary and has been previously seen as a spontaneous feature that can destabilize some stable systems. Therefore, differential equations with time-delayed influences have attracted much interest in current history since they have been demonstrated to be helpful resources [1][2][3] . Nonetheless, time-delay has been used lately in a number of research papers as a controlling factor that can be used to develop the regulator effectiveness. Actually, numerous physical processes have been represented by employing nonlinear ordinary differential equations in the previous decades. The DO is just one of these equations that have attracted great attention because of its conventional uses in various areas. Several researchers are motivated by this nonlinear differential equation because it reflects the normalization approach in our real world, based on its unique oscillating and chaotic environment [4][5][6] . The DO arises when a body moves under the influence of nonlinear spring strength, linear stickiness dampening, and periodic powers 7 . It could be used to establish mechanical experiences and significant changes by periodical exterior strength. An effort to explore several DO applications and propose an alternate mathematical model to approximate them was made 8 . In this approach, the power series was implemented as the beginning function through the integration brought about by a one-step interval. With regard to calculations of circuit analysis, the DO is supplementary and particularly difficult to deal with. Therefore, Chen et al. 9 scrutinized the circuits from the standpoint of vibrating sciences and energetics available to support the investigation of their properties. As a result of this work, users will be able to comprehend the mechanics and properties of DO in a much more straightforward way. The influence of an increased simulation on significant stationary-state formulations and bi-stability was investigated 10 in a delayed DO. The fundamental equation of the slow dynamical oscillator and the related flow was obtained using the honest partitioning of motion approach and the multiple scales method 11 . Numerical computations have validated the quantitative outcomes of this study. This method was used in many works [12][13][14][15] to attain the estimated formulas of various vibrational models. The stability of the parametric DO was investigated 16 . The decelerated factor and the cubic stiffening factor were shown to have a destabilizing effect. Fundamental and parametric harmonics have stabilizing effects. The cubic-quintic nonlinear term and the externally stimulated impulse were explored in a combined Rayleigh-Van der Pol-Duffing oscillator 17 . To provide an approximate solution, the Poincare'-Lindstedt methodology was modified 11 . A comparison of the estimated approach and RK4 has revealed significant compatibility.
Differential equations are widely utilized in computational chemistry in addition to industrial biomechanics in developmental stages. Resolutions towards such equations have been required to assess and enhance reactions in a mathematical conception. Computing an accurate solution is challenged in the majority of situations. To understand such difficulties, a variety of mathematical approaches has been modeled. Therefore, simplifications and computational methodologies are widely employed since most equations may not have an accurate solution. Subsequently, many classical perturbation methods were recognized. Unfortunately, all of these methods necessitate the existence of a small parameter (SP) in the handheld equation. The establishment of a synthetic SP is frequently motivated by the lack of a meaningful physical application of a zone to be too limited. In a few words, the SP may be applied in a variety of ways. So, the supplementary critical question arises. Is it imperative to attain information immediately? It may be done by using SP or using a synthetic one. In this regard, it is preferable to talk properly about "methodologies specialized in parameter improvement instead of merely the SP. Really, there is no distinction between real and artificial SP from this viewpoint. Nevertheless, the expression manufactured SP should be used again as is customary. It is worth mentioning that the presence of SP has been suggested in a diversity of mathematical fields. The Chinese Mathematician Prof. He [18][19][20] created the HPM in response to the constraints of the classical perturbation methods. As an essential concept, the SP has no significance on the HPM. In accordance with the Homotopy equation, the given problem may be constructed into two linear and nonlinear parts and divided by the parameter ρ ∈ [0, 1] . It has been demonstrated that it can be used to analyze highly nonlinear vibrational issues and create a multiplicity of real-world engineering challenges. The fundamental equations of exceedingly nonlinear excitations of orthotropic structural components are deduced using the Von Karman huge defection hypothesis and the Galerkin method to achieve a precise and effective alternative approach to the nonlinear vibrational problem of the membrane 21 . Using HPM, an analytical solution to the problem of free vibration of rectangular plates with convoluted border circumstances in a plane was suggested 22 . The HPM was used to introduce a solution to a time-delayed differential equation and was followed by numerical illustrations 23 . The outcomes demonstrated that the suggested strategy is both successful and straightforward to implement. For an active vibration control of a beam on an elastic base, the HPM was used 24 . The methodology was used to load an axially loaded elastic basis with one fixed end, while the other was simply maintained. In comparison to other interpolating techniques, it has been shown that HPM provides a highly accurate solution, in association with the other techniques, for a specific problem.
The concept of control was widely incorporated into the modeling of autonomous engines or physiological systems, even though the feedback was required to maintain a stable state. But a lot of this feedback required a short amount of time to analyze the data. It was common practice to describe this occurrence by formulating a delay differential equation (DDE), in which the evolution of a dependent variable at a time t depends on its value at a time (t − τ ) . Unfortunately, it was a very difficult mathematical problem to solve a DDE. In the last years, there has been a significant increase in processing power, which has reignited interest in DDEs 25 26 . A class of nonlinear differential equations was studied using this novel formula. The Adomian solution utilizing this new formula converges more quickly, according to a numerical experiment. Almost all scientific disciplines, including mathematics, physics, engineering, biology, neuroscience, physiology, economics, and many more, are connected by the topic of dynamic systems with time delay, which was an active area of research. Together, contributions from theoretical and experimental groups were highlighted, stressing a wide range of applications 27 . There were reports on systematic bifurcation studies of mode-locked laser systems, higher-order locking phenomena between delay and laser oscillation period, and management of cavity solitons as light spots in spatially extended systems. Both experimental and analytical methods were used to examine the nonlinear recovery of reverse-biased waveguide absorbers based on quantum dots 28 . The recovery dynamics were demonstrated to consist of a quick initial layer and comparatively moderate decay. The carrier capture, escape, and Pauli blocking processes were taken into consideration in a rate equation model of a quantum dot semiconductor optical amplifier 29 . It examined potential distinctions between Auger processes that are prominent for recovery and phonon-assisted processes. Experimental, numerical, and analytical research were done to determine how a nonlinear optical oscillator will react to delayed broadband bandpass filtering feedback 30 . The examination of the DO, in light of the previous advantages, has potential requests in a wide range of physical and technological phenomena. Therefore, continuing employment studies on the stability analysis of the parametric DO have an extensive variety of applications in science and engineering. A second-order ordinary differential equation with a cubic-quintic nonlinearity is represented in this equation. To establish a better approximate solution, the HPM is joined with the Laplace transform, and the idea of expanding frequency will be utilized. A comparison with RK4 is made to verify the prototypical association. Consequently, the existing work mainly discussed DO and studied its excitation with location and velocity delays in the following way: where all of the used constants in Eq. (1) are previously provided.
To crystallize the implementation of our inquiry, the rest of the paper will be outlined along with the following Sections: an improved approximate solution to the considered problem using the improved HPM is explained in Section "Improved HPM". Section "Discussion of the HPM outcomes" is devoted to the discussion of the previous outcomes. The mathematical procedure of the nonlinear stability analysis is introduced, and the primary resonance is established in Section "Mathematical Formulation of nonlinear analysis". The frequency response on numerical validations is discussed in depth in Section "The frequency-response curves and numerical validations". Additionally, its subsections examine the controlled system with/without the time delay as well as the uncontrolled system. Finally, the main findings of the current study are summarized in "Conclusion" Section.
It is easier to imagine the initial circumstances as follows: The HPM technology is constructed mostly on the underlying Homotopy equation:

Improved HPM
The standard method could be changed to respond accurately over time. Such a skillfully adapted HPM could very well reveal the previously underappreciated importance of the time delay parameter. Nevertheless, it also addresses the problem by using an exponential delay to offer a corresponding characteristic across the decay constriction. One can get a constrained solution over time, and that is the topic of this essay. To explore the effects of the delay parameter, we might evaluate Homotopy Eq. (3) using the new expansion that departs from the traditional one. In conclusion, the time-dependent amount can indeed be increased in the way described 31 : where y j (j = 0, 1, 2) are the stages of y . The recognized natural frequency can be enlarged as follows using the methodology of the previously understood investigation 16 : The developing formula that is provided in Eq. (4) is inserted into the Homotopy equation, in light of the starting conditions as given in Eq. (2), to generate an improved approximate solution of Eq. (3). As an outcome, the analytic mathematical formulation of the zero-order problem is specifically defined as follows: The first problem of the Homotopy Eq. (3) can be described as: www.nature.com/scientificreports/ The regularly homogeneous improvement requires the removal of the secular associations when Eq. (6) is involved in Eq. (7). In general, the coefficients of the circular functions sin ̟ t and cos ̟ t reflect the behavior of the fundamental regular situations. Therefore, the uniform valid approximate solution needs the following restrictions: and The solvability criteria are specified by Eqs. (8) and (9), respectively. Our aim, away from the presence of secular terms, is to obtain a consistently reasonable solution for the firstorder stage y 1 . At this point, the first-order Eq. (7) then becomes Taking Laplace transforms for both sides of Eq. (10), we attain the following solution One more, solving Eqs. (8) and (9), we achieve and To this order of approximation, Eqs. (12) and (13) are combined, keeping in mind ω 1 = ̟ 2 − ω 2 , to construct the characteristic equation of the expanded nonlinear frequency as: To this end, the periodic estimated solution of the prepared principle of movement in Eq. (1) can be represented as follows:

Discussion of the HPM Outcomes
This section is devoted to graphing and analyzing the achieved outcomes by using the improved HPM. To confirm the approach of our procedure, it is desirable to display both the analytical and numerical solutions. A prototype specified selection system that has the following characteristics is used: Therefore, the analytic solution (AS) is graphed in Fig. 1 when the decay parameter has the following different amounts τ (= 0.5, 0.8, 1.1) . The subsequent figure is calculated when one of the real roots of Eq. (14) is considered. With the aid of the Mathematica Software, we choose the real reasonable value of the nonlinear expanded frequency as ̟ = 6.84115 . As shown, the curves of Fig. 1 have an oscillatory decay form till the end of the examined time interval which is consistent with Eq. (15). When τ increases, the oscillations decay increases, and the curves have a stationary performance at the end of the time interval. The corresponding phase plane diagrams of these curves are drawn in Fig. 2. It is observed that the curves start from one point, which is the initial point, and have spiral forms until the end at another stable point. This signifies that the obtained solutions are stable.
The numerical solution (NS) of Eq. (1) is obtained by employing the RK4 and is compared with the AS when the previous data are used at τ = 0.5 as shown in the curves of Fig. 3, while the phase plane curves are drawn in Fig. 4. The difference between them is due to the forms of Eqs. (1) and (15). The curves of Fig. 5 reveal the (15) y = lim ρ→1 e −ρtτ (y 0 (t) + ρy 1 (t)).

Mathematical formulation of nonlinear analysis
Throughout this Section, forced and damped vibrations are considered by means of the multiple scale homotopy perturbation method. The solutions will be verified with the numerical imitations. Presume the following hampered structure with exterior excitement.    Fig. 6. This figure shows how the A/D process can be executed sequentially, where the macro-fiber composite (MFC) sensor measures the oscillations of the system. It follows that the measured signal is fed to a digital signal processor (DSP) that acts as a controller via an A/D. The previous control algorithm performs the required mathematical processing for the acquired signal. The processed signal is then fed back via a power amplifier to the MFC actuator which, in turn, improves the vibrational behavior of the system. Therefore, to compare the efficiencies of the proposed controllers and to report the optimum control factors, a nonlinear examination of the whole structure (i.e., Eqs. (16) and (17) is established. However, to choose the optimal control factors ( γ , , and τ ) that can minimize the oscillation amplitude of the assumed scheme, a scientific examination of the impact of these parameters on the scheme of dynamical performance must be presented by evaluating the scheme of dynamical Eq. (1). Whereas the hampering and exterior excitement with the control parameters are rearranged so that they can balance the impact of nonlinearity on the maximum order. This choice is particularly appropriate for the basic resonances of oscillations. (16) y + ω 2 y + 2µẏ 2 + αy 3 + βy 5 = F cos �t + F c ,  where ρ is the embedded artificial homotopy parameter, and L(y) and N(y) are the linear and nonlinear portions of the assumed differential equation, correspondingly. They are identified as: and Consequently, the homotopy formula can be constructed as: Without sacrificing generality, two-time scales could be contemplated. Characteristically, the shift may be extended as follows: (27) y 0 (T 0 , T 1 ) = A(T 1 )e iωT 0 + cc, www.nature.com/scientificreports/ Let σ denote a detuning parameter that expresses the nearness of the excitement frequency to the scheme natural frequency ω such as (� = ω + ρσ ) . To attain an unbroken effective expression, the secular terms should be excluded. The removal of these parts needs the termination of the quantities of the functions e ±iωT 0 . Accordingly, one derives the subsequent solvability situation as follows Appropriately, the solution of Eq. (28) must be formulated as: Multiplying Eq. (29) by ρ , and then considering the tendency to unity yield. Equation (31) denotes a first order nonlinear differential equation with a complex coefficient; the solution of the function A(t) can be stated in the following polar form: where ψ(t) and a(t) are two real functions concerning time.
The above functions correspond to both the improved phase-angle and the vibration amplitude of the scheme, correspondingly. Inserting Eq. (32) in Eq. (31) and dividing the real and the imaginary elements, one finds the following amplitude-phase adjustment equations: and where φ = σ t − ψ.
At steady-state vibrations, one gets ȧ =φ = 0 27 ; therefore, Eqs. (33) and (34) become and The combination of Eqs. (35) and (36) yields In the stability diagram, the amplitude vibration will be graphed against any factor of the scheme like σ , F, γ , or τ . Using Mathematica software, the following diagrams will be demonstrated soon in the following section. As opposed to that, the linearized stability can be investigated near the fixed points (equilibrium points). Therefore, one can use the following consistent assumptions 33,34 : Inserting Eq. (38) into Eqs. (33) and (34), and retaining just the linear elements, a little departure from a 11 and φ 11 is obtained. The following standard expression is the outcome of this approach Established on the last square matrix (Jacobian matrix), the stability of Eq. (39) can be investigated by examining the Jacobian matrix eigenvalues. (28) F cos φ.
(39) Therefore, the following can be expressed as the required and adequate requirements of stability:

The frequency-response curves and numerical validations
This section investigates the operational characteristics of the model under consideration by resolving the system frequency-response equations at the fundamental resonance situation 35 . The assumed outcomes are attained based on the selection of the dimensionless system data as: ω = 4.5, µ = 0.05, α = 0.25, β = 0.5,F = 0.1, γ = 0, = 0, and τ = 0 . Furthermore, the numerical approach to the used parameters yields different figures explaining the structure of the fundamental equation Eq. (1) with the help of the Mathematica Software, at various amounts of the equivalent parameters σ , F, γ , , and τ . The stationary-situation of the vibration amplitude is numerically attained. The next subsections are structured so that subsection "Unrestrained structure" can examine the oscillatory behaviors of the unrestrained scheme, and subsection "The controlled system without time-delay" compares the effectiveness of the two suggested controls at reducing vibrations while ignoring loop delay. The final part "The controlled system with time-delay" investigates how the two controls behave dynamically when the loop delay is reflected.
Unrestrained Structure. The nearness of the excitement frequency to the structure natural frequency ω is formulated using a detuning factor σ . Consequently, one may employ σ as a controller factor to examine the system frequency-response curve in the main resonance situation. Figure 7a examines the frequency-response curves of the uncontrolled system at various stages of the excitement strength F . The plotted curves in this figure illustrate that the vibrations of the system monotonically increase with the increase of this force. Furthermore, the scheme may be excited with the amplitude of the considerable oscillation, which depends on both the excitement frequency and the excitement amplitude F . On the other hand, curves of Fig. 7b represent the numerical simulation of the oscillations time histories which corresponds to the curves in Fig. 7a at σ = 0 (i.e., . (44) (45) The controlled system without time-delay. As shown in the following discussion, the impacts of the control factors γ and on the scheme response curves are plotted in Fig. 8. Figure 8a demonstrates the influence of the situation gain γ on the frequency-response curves of the system. The linear controller hypothesis makes it clear that the position of the feedback controller is what causes changes to the linear natural frequency system. Accordingly, the frequency-response curve shifts as the positional gains grow. In this subsection, the effect of the two controls gains on the oscillatory behaviors of the scheme is examined when the loop delay vanishes. It is clear from Fig. 8a that, based on γ values, the linear-position controller is accountable for any variation of the scheme frequency response curve. Consequently, it is preferable to ignore the amplitude high oscillation of the system with the variation of its normal frequency, away from the excitement frequency, by using the linear position. The numerical reproductions, for the controlled scheme temporal vibrations according to the plotted curves in Fig. 8a at σ = 0 (i.e., � = ω ), are graphed in Fig. 8b. The numerical solutions of the original Eq. (1) are in high agreement with the behaviors of the curves included in Fig. 8a. Moreover, the oscillations of Fig. 8a have some convergent beams and a stable manner to the end of the examined time interval.    Fig. 9a show the influence of decreasing the linear-velocity feedback gain on the scheme frequency-response curve, where the decrease of improves the structure of the linear damping parameter, which eventually decreases the scheme vibration amplitudes. It is observed that with the increase of , the frequency response curve increases, where three pick points are observed with the used values of . The nonlinear temporal vibration of the beam system, when decreasing the velocity feedback gain from = −0.1 to = −1 according to Fig. 9a, is numerically proven in Fig. 9b. The comparison between the behaviors of both figures reveals good consistency between their behaviors.
The controlled system with time-delay. To discover the impact of the controller loop-delay τ on the stability and the vibration amplitude of the system, curves of Fig. 10a are plotted. The position gain γ is employed as a controller factor against the time delay τ whenever σ = = 0 to describe the stable and unstable zones in (γ − τ ) plane. This figure demonstrates that the unstable zone of the solution is a periodic function of time-delay either for negative or positive position gain. Consequently, it is likely to choose the loop delay in a particular approach for a positive or negative position gain to ensure the stability of the system. Figure 10b,c display the system frequency-response curve at various amounts of γ and τ which have been selected according to Fig. 10a. Figure 10b  www.nature.com/scientificreports/ τ = 0.3 and τ = 1.5, respectively. Consequently, one can improve the vibration decreasing efficiency of the control employing the loop-delays as a new controller factor, where the optimal amounts of the loop-delay τ that can improve the control efficiency must be chosen within the middle of the stable solution zone. At τ (= 0.3, 1) , the system frequency-response curve has two unstable regions in (γ − τ ) plane at γ (= −1, 1) , where the structure will achieve unstable motion for any excitation frequency. To confirm the precision of the achieved stability chart in Fig. 10a, a numerical imitation for the deemed scheme has been graphed in Fig. 10d. The comparison of the curves included in Fig. 10b with Fig. 10d displays that the structure replies with bounded oscillation amplitude as τ (= 1, 1.35) because these points are chosen inside the stable zone of the solution. As opposed to that, the regulated scheme replies with boundless oscillation amplitude at τ = 0.5 as this value has been chosen inside the unstable solution zone. The comparison between Figs. 10c,d has the same direction as before with Figs. 10b,d.   Fig. 11a illustrate the stable and unstable solution areas in the plane ( − τ ) when σ = γ = 0 . The graphs demonstrate that these zones are periodic functions of either positive or negative velocity gain in terms of time-delay. The structure frequency-response curve at various points in ( − τ ) plane is depicted in Fig. 11b-e. Furthermore, the curves of Fig. 11b describe the frequency-response curve of the scheme and are computed when the velocity feedback gain = −1 at τ (= 0.1, 0.7, 1.3, 1.4) to show the influence of the loop-delay on the scheme, where the best time-delay at speed feedback gain is found at τ = 1.4.
The three values of delay τ (= 0.1, 1.3, 1.4) have the same amplitude at σ = 0 , which is numerically in an agreement with Fig. 11c. The controlled scheme replies with a boundless oscillation amplitude at τ = 0.7 and it is of the same consistency with the corresponding curve behavior of Fig. 11c at the same value of τ . When τ (= 0.1, 0.7, 1.3, 1.4) the scheme frequency-response curve of speed feedback gain = 0.2 is demonstrated in Fig. 11d, where the best time-delay for vibration destruction is at τ = 0.7 . Based on Fig. 11a, choosing the loopdelay to be at the center of the stable solution zone in ( − τ ) plane improves the effectiveness of the controller. Moreover, Fig. 11d has two stable and two unstable zones in the plane ( − τ ) at (= 0.2, −0.2) and τ (= 0.7, 1.3) , where the structure will operate unstable motion for any excitation frequency. Figure 11d shows high consistency with Fig. 11e.
It is likely to prevent the structure high oscillation amplitude with the variation of the linear-position and the linear-velocity feedback control. As shown in Fig. 12a, the linear-position with linear-velocity feedback control is more effective with time delay in reducing the structure nonlinear vibrations. It is described that the optimal time-delays amounts can increase the structure strength against instability based on the obtained analytical and numerical investigations in Fig. 12b. Furthermore, the system vibrations should be chosen so that they maximize the linear damping function µ Eq = µ + (γ ω) sin(ωτ ) − cos(ωτ ) to increase the linear (position-velocity) controllers efficiency in the suppression case. Additionally, it is found that the change in time delay suppresses the system vibration. Based on the curves of Fig. 12a, we can verify that the curves in Fig. 12b represent the optimal control strategy for reducing the extreme oscillation amplitude of a parametrically excited structure as well as feedback gains (i.e., γ = 2, = 0.8).

Conclusions
The subject of this paper is the time-delayed controllers of a damped nonlinear stimulated DO. Since timedelayed techniques have lately attracted a lot of research attention, the investigation subject is very contemporary. Time delays of position and velocity are used to reduce the nonlinear oscillation of the existing model. A modified HPM is employed to obtain a significantly more precise approximate solution [36][37][38] . For various amounts of the used factors, the temporal difference of this solution is graphed. The analytical approximate solution has been compared with the numerical solution, which has been obtained by employing the RK4 and they have proved consistent. The results have demonstrated that the improved HPM is suitable for a range of damped nonlinear oscillators because it reduces solution error while boosting validation diversity. Additionally, it offers a viable model that addresses a range of nonlinear issues. Gained has been the nonlinear algebraic equation that controls the steady-state oscillation amplitude. Investigations have been done on the effectiveness of the suggested controller, the influence of time delays, control gains, and feedback gains. The actual results have demonstrated how the sum of the time delays in the control loop and the product of the control and feedback gains affect the controller performance. Analytical and numerical analyses have shown that optimizing the suggested controller could entirely cut down on system vibrations for particular values of the control and feedback signal. It is possible to obtain the nonlinear algebraic equation governing the steady-state oscillation amplitude. When there are time delays as well as control and feedback gains, the controller performs well. Though the method adopted www.nature.com/scientificreports/ in this study is traditional, it has yielded novel results. The findings have shown that the sum of time delays in the control loop as well as the product of control and feedback gains have an effect on the controller effectiveness. Theoretical and numerical models have demonstrated that the suggested control may successfully lessen structure vibrations at suitable levels of the controller and feedback signal gains.